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Abstract 

We consider a (2 + l)-dimensional field theory, assumed to be holographically dual to the extremal 
Reissner-Nordstrom AdS4 black hole background, and calculate the retarded correlators of charge 
(vector) current and energy-momentum (tensor) operators at finite momentum and frequency. We 
show that, similar to what was observed previously for the correlators of scalar and spinor opera- 
tors, these correlators exhibit emergent scaling behavior at low frequency. We numerically compute 
the electromagnetic and gravitational quasinormal frequencies (in the shear channel) of the extremal 
Reissner-Nordstrom AdS4 black hole corresponding to the spectrum of poles in the retarded correla- 
tors. The picture that emerges is quite simple: there is a branch cut along the negative imaginary 
frequency axis, and a series of isolated poles corresponding to damped excitations. All of these poles 
are always in the lower half complex frequency plane, indicating stability. We show that this analytic 
structure can be understood as the proper limit of finite temperature results as T is taken to zero 
holding the chemical potential \i fixed. 



1 Introduction 



The AdS/CFT correspondence is emerging as a promising approach for analyzing strongly-coupled 
condensed matter systems. 1 Although the systems accessible through the duality need a lot of refine- 
ment in order to resemble actual condensed matter systems, (more than) a decade-long experience 
with the correspondence has taught us that it is extremely useful in extracting some universal features 
of strongly-interacting systems which would be hard, or even impossible, to obtain using field theoretic 
methods. It is in this sense that applying the tools of the correspondence to study strongly-interacting 
condensed matter systems seems to be promising. Recent examples of such applications include mod- 
eling superconductivity and superfluidity [2, 3, 4, 5, 6], engineering systems with Schrodinger and 
Lifshitz symmetries [7, 8, 9, 10, 11, 12, 13, 14], or modeling a system which can exhibit non-Fermi 
liquid type behavior [15, 16, 17, 18, 19]. See [20, 21, 22] for reviews. 

Using the correspondence to study quantum critical regions of condensed matter systems is another 
example of such applications [23, 24, 25]. At the quantum critical point, the system is often a strongly- 
coupled conformal field theory, and it is usually the case that a proper field theoretic understanding 
is not known. The authors of [17] considered a simple gravitational background whose dual field 
theory exhibits a variety of emergent quantum critical behaviors. They considered the background 
of extremal Reissner-Nordstrom AdS^+i black hole whose dual field theory is supposed to be a d- 
dimensional strongly-coupled field theory which is at zero temperature and finite U(l) charge density. 
Considering probe scalars as well as spinors in the geometry, it was shown that the retarded correlator 
of the dual operators show emergent quantum criticality at low frequency. The existence of such 
behaviors was attributed to the fact that the near horizon region of the background geometry (which 
encodes the IR physics of the boundary theory) is AdS2 x B 2 , and the assumption that there exists 
some sort of IR CFT dual to the AdS2 region. It was shown that the behavior of the retarded Green 
function Gr(ui, k) at low frequency is dictated by Gr(uj = 0, k) as well as the conformal dimension of 
some operators in the IR CFT. More specifically, it was shown that Gr(uj, k) can in general exhibit a 
scaling behavior for the spectral density, a log-periodic behavior, or in the case of fermionic operators 
indicate the existence of Fermi-like surfaces with quasi-particle excitations of non- Fermi liquid type. 

In [26] we examined the role played by the IR CFT in constraining the low frequency behavior of 
retarded correlators of charge vector and energy momentum tensor operators of the boundary the- 
ory. The main focus though was on extracting some transport coefficients, such as conductivity and 
viscosity, of the boundary theory, and investigating the relationship between the IR CFT and the uni- 
versality of such coefficients. As such, the analyses were limited to the case of zero spatial momentum. 
The boundary theory considered was a (2 + l)-dimensional field theory at zero temperature and finite 
charge density dual to an extremal Reissner-Nordstrom AdS4 black hole. Although temperature is 
zero, there is a scale in the problem set by the chemical potential for the charge density and moreover 
the entropy density is finite. 2 It was shown in [26] that, for example, the ratio of shear viscosity to 
the entropy density is l/4ir saturating the KSS bound [27]. There are a variety of reasons why such 
analyses at zero-temperature are interesting. First, calculating transport coefficients involves taking 
the limit of small frequency, which at zero temperature becomes subtle and extra care is needed to deal 
with this subtlety. Secondly, they generalize the universality arguments, given for finite temperature 

1 See [1] for a review of the AdS/CFT correspondence. 

2 This is presumably a "large N" effect. Analysis of this has recently appeared in [19]. 
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backgrounds [29, 28, 30], to backgrounds of extremal black holes. Third, the system studied may be 
in the universality class of some quantum critical points, hence understanding the universal behavior 
of the transport coefficients of this system may shed light on the physics of quantum critical points. 

In this paper, we analyze the retarded correlators of the charge vector current and energy-momentum 
(tensor) operators of the aforementioned (2 + l)-dimensional boundary theory for non-vanishing (spa- 
tial) momentum. Although we perform the analyses in the (2 + l)-dimensional boundary theory, 
they can easily be carried out for higher dimensional boundary theories dual to extremal Reissner- 
Nordstrom AdS^+i black hole. Our focus here is on the retarded correlators of those operators which 
are in the so-called shear channel of the boundary theory. We find that for generic momentum, similar 
to the cases of scalar and spinor operators studied in [17], the retarded correlators of such vector and 
tensor operators exhibit emergent scaling behavior at low frequency. These operators give rise in the 
IR to two sets of operators with different conformal dimensions. The low frequency scaling behaviors 
of the boundary theory retarded correlators, especially their scaling exponents, depend only on one of 
the two IR CFT operators. By taking linear combinations of the dual bulk modes, we construct the 
gauge invariant modes appropriate for extracting the quasinormal modes. 

We study the small frequency regime analytically for finite momentum, and we explore the ana- 
lytic structure of the shear channel retarded Green functions numerically for generic frequency and 
momentum. We present evidence for the existence of a branch cut in these Green functions at the 
origin, as well as a series of metastable modes corresponding to isolated poles in the lower half complex 
frequency plane. By turning on temperature we observe how the branch cut dissolves into a series 
of poles on the negative imaginary axis and determine the corresponding dispersion constant of the 
leading pole, in agreement with previous results. 

The organization of the paper is as follows. In section 2, we review the Reissner- Nordstrom AdS4 
black hole background, its dual boundary field theory, near horizon AdS2 xR 2 geometry and the dual IR 
CFT. In section 3 we present the Linearized Einstein-Maxwell equations for the bulk modes dual to the 
(vector and tensor) operators in the shear channel of the boundary theory. Taking linear combinations 
of the modes, we construct the gauge- invariant combinations in terms of which the linearized Einstein- 
Maxwell equations reduce to a set of two coupled second-order differential equations. By introducing 
"master variables" we then decouple these equations. In section 4 we solve the decoupled equations 
by matching the solutions in the outer region to the solutions in the inner AdS2 region. Having 
determined the solutions, we calculate the retarded correlators of the charge vector current and energy- 
momentum tensor operators of the boundary theory. For generic momentum, we then extract the low 
frequency emergent scaling behavior of their spectral functions. In section 5 we argue that the (shear- 
type) electromagnetic and gravitational perturbations of the extremal Reissner-Nordstrom AdS4 do 
not cause instability by showing that the associated quasinormal frequencies are all located in the 
lower half of the complex frequency plane. We then numerically compute the spectrum of the above- 
mentioned quasinormal frequencies and compare the numerical results with the analogous quasinormal 
frequencies of the non-extremal Reissner-Nordstrom AdS4 black hole background. In Appendix A we 
have summarized the finite-temperature equations we occasionally used in the bulk of the paper in 
order to extract the shear-type quasinormal frequencies of the non-extremal Reissner-Nordstrom AdS4 
black hole. 
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2 The Background and the Boundary Field Theory 



Consider the Einstein-Maxwell action in 3+1 spacetime dimensions with a negative cosmological con- 
stant A = -3/L 2 

S = ^2 [ d 4 x^ (R-2A- L 2 F^) , (2.1) 

where L is the curvature radius of AdS4. The background we consider is the Reissner-Nordstrom AdS4 
black hole 



ds 



2 



A 



9iiv dx»dx v = ^ ( - f(r)dt 2 + dx 2 + dy 2 ) + ^J^dr 2 , (2.2) 
^(l- r f)dt, (2.3) 



which is a solution to the Einstein-Maxwell equations obtained from (2.1), where 

/<r) = l- if „ = ^>. (2.4) 

Here, tq is the horizon radius and is given by the largest real root of /(ro) = 0. Also, note that 
M = 1 + Q 2 . The temperature of the black hole takes the form 

while its entropy, charge and energy densities are given by 

respectively [31]. The background is extremal when Q 2 = 3, for which the temperature vanishes but 
the entropy density remains finite. The background is invariant under flipping the sign of At. We 
choose fx to be positive resulting in Q to be positive. With this convention, Q = y/3 at extremality. 

The background is dual to a (2 + l)-dimensional strongly-coupled field theory at finite temperature 
T and finite charge density p. The entropy and energy densities of the dual theory are given by s and 
e in (2.6), respectively. Also, the asymptotic value of the bulk gauge field A t (co) = il is interpreted in 
the dual theory as the chemical potential for the (electric) charge density. Little is known about the 
details of the dual theory from field theory perspectives. On the other hand, using holography a lot 
has been learned (especially thermodynamical properties) regarding its strong-coupling behavior; see 
[31] and its citations. 

In this paper we work in the extremal limit, where the dual theory is at zero temperature but finite 
charge density. We will refer to this dual theory as the boundary field theory. Although the black hole 
temperature vanishes at extremality, its horizon area remains finite whose dual interpretation is that 
the boundary theory has a finite entropy density at zero temperature. 

In the extremal limit, /(r) in the background metric (2.2) takes the form 

«r)-l-«(2)* + »(2) 4 , (2.7) 
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which has a double zero at the horizon, and can be approximated near that region (to the leading 
order in r — tq) by 

/(r)^-|(r-r ) 2 . (2.8) 

The near horizon geometry is AdS2 x M 2 . To see the emergence of this geometry, first change the 
radial coordinate r to r] defined by 

r~r = ^. (2.9) 

There is then a scaling limit [17] in which 

ds 2 = ^ ( - dt 2 + drj 2 ) + T f 2 (dx 2 + dy 2 ) , A = ^dt. (2.10) 

The curvature radius of the AdS2 is L2 = L/y/Q. The radial coordinate is interpreted holographically 
as the renormalization scale of the dual field theory, and the near horizon region corresponds to the 
IR limit. This implies that the AdS2 x M 2 geometry encodes the IR physics (00 — > 0) of the boundary 
theory. 

On general grounds of holography, one expects the gravity on the AdS2 space to be dual to a 
CFTi. This led the authors of [17] to suggest that the (2 + l)-dimensional boundary field theory 
(which is dual to the extremal charged AdS4 black hole) flows in the IR to a fixed point described 
by a CFTi. Following [17] we will refer to this CFTi as the IR CFT. The details of the AdS 2 /CFTi 
correspondence and how exactly the mapping works are poorly understood. In particular, it is not 
clear whether the CFTi of the correspondence represents a conformal quantum mechanics or a chiral 
sector of a (l+l)-dimensional CFT. Nevertheless, not knowing the details of the IR CFT and just 
assuming its existence, the authors of [17] were able to show that the low frequency behavior of some 
observables in the boundary CFT is encoded in the IR CFT. In what follows, we consider the charge 
current and the energy-momentum tensor operators of the boundary field theory and elucidate the 
role of the IR CFT in determining the low frequency behavior of the retarded Green functions of the 
two operators. 



3 Gauge Field and Metric Fluctuations 

The charge current operator in the boundary field theory is dual to the fluctuations of the bulk gauge 
field in the transverse direction while the energy-momentum tensor operator is dual to the fluctuations 
of the background metric. The calculation of the corresponding retarded Green functions starts with 
solving the linearized Einstein-Maxwell equations for these fluctuations. 

3.1 Einstein-Maxwell Equations 

To obtain the linearized Einstein-Maxwell equations, we first define 

9iiv = 9iiv + h^ u , A^ = A^ + a^, (3-1) 
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where and represent, respectively, the (extremal) metric and the gauge field of the background, 
and h^u and a M are the fluctuations. We choose the so-called radial gauge 



where u = {t,x,y,r}. We proceed by Fourier transforming the fluctuations 



(3.2) 



(3.3) 



where, without loss of generality, we used the rotation invariance in the (x, y) plane to set k y = and 
defined k x = k. The fluctuations split into decoupled groups depending on whether they are even or 
odd with respect to parity, y — > —y. Accordingly, h ty , h xy , a y have odd parity while ha, h tx , h xx , h yy , 
at, a x all have even parity. In this paper, we study the odd parity modes which translates into the 
shear and (charge) diffusion modes of the boundary theory. The analysis of the even parity modes 
(sound modes) will be given elsewhere. 

It is more convenient to raise the indices in ht y and h xy (by the background metric g^u) and work 
with h y t, h x y . It is also convenient to define the dimensionless quantities 



u 



to 



k 

3 



(3.4) 



In this notation, the linearized Einstein-Maxwell equations for the odd parity modes h* y , h x y , a y then 
read 

r 12 i 

f(u) u 4 h y 'l(u) + 4u 3 h y, t (u) + — a' y (u) - 3q [mh x y (u) + qh y t (u)] = 0, (3.5) 

L /i J 

f(u) [u A f(u)h^(u) + [u 4 f (u) + 4n 3 /(n)] h x, y (u)\ + 3to [voh x y (u) + qh y t (u)} = 0, (3.6) 
f{u) [« 4 /(«)a» + u 2 [u 2 f(u) + 2uf{u)] a' y (u) + [m 2 h y ' t (u)] +3[vo 2 - /(n)q 2 ] a y {u) = 0, (3.7) 

where f(u) is given in (2.7). There is also a constraint which comes from the yn-component of the 
linearized Einstein equations given by 



u A vo h y ' t (u) + u 4 f(u)q h x ' y {u) + —to ay(u) = 0. 



(3.8) 



As we will see later in more detail, this constraint, evaluated asymptotically, encodes a Ward identity 
in the dual boundary theory, as one should expect. 



3.2 Gauge Invariant Modes 

Although working in the radial gauge (3.2) is convenient, it does not completely fix the gauge freedom 
of a M and h^ u . In the present case, it is easy to identify the gauge invariant combinations of fields. 

For the metric and gauge field fluctuations of our background, there is a U(l) gauge transformation 
An — > An = Ap + V^X) according to which f^v = d^a-v — d^a^ is invariant. There are also gauge 
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transformations associated with diffeomorphism under an arbitrary vector field £, implying that the 
fluctuations {/j MJ ,,a p } and {h^ v ,a p } are equivalent if they are related as [32] 



h 



III/ 



(3.9) 
(3.10) 



where V is any torsionless connection and V is the connection associated with g pu . It is important to 
notice that these gauge transformations are not changes of the coordinates. In particular, they affect 
only the fields and a M , while all the other tensors, vectors, etc. remain invariant. 

In this paper, the fluctuations of interest are h y t, h x y and a y . While a y is invariant under the residual 
gauge transformations, h y t and h x y transform according to h y t — > h y t+i^ivo^ y and h x y — > h x y — ifi<\^ y , 
respectively. Thus 



X(u) = qh y t (u) + voh x y (u), 
Y{u) = a y {u), 



(3.11) 
(3.12) 



are invariant combinations. 



Note that at q = the symmetries are enhanced, and h x y {u) is gauge invariant. For this special 
case, as can be seen in equation (3.6), h x y (u) decouples from the rest, while equations (3.7) and (3.8) 
together give a decoupled equation for a y (u). The decoupled equations read 



3t» 5 



u 4 f(u)h^(u) + [u 4 f{u) + 4u 3 / («)] h x ' y (u) + J^h%(u) = 0, 



12 



V(u)a» + [u 2 f'(u) + 2u/(u)] a» + ^ " ^ 



0. 



(3.13) 
(3.14) 



Equations (3.13) and (3.13) were studied in [26] and the retarded Green functions of the dual operators, 
T xy and J y , at q = and small vo were calculated. For the most part in this paper, we consider q to 
be finite and non-zero. 

Even at finite q and vo, the Einstein-Maxwell shear channel equations can be decoupled in terms 
of a pair of "master fields" [33]. This is remarkable, since in the boundary field theory it corresponds 
to diagonalizing the entire renormalization group flow. The master fields are 



vo 2 — /(n)q 2 u [vo 2 — / (u)q 2 

and they satisfy the decoupled equations 

3 / o „, s o\ 6 



+ u 1± 



Y(u). 



(3.15) 



[u 2 f( U )$> ± (u)]' + uf(u) + ^ (vo 2 - f(uW) -- 3 (i± yrr^) 



$ ± (u) = 0. (3.16) 



As we will discuss in the next section, it is most convenient to organize numerical calculations in terms 
of these modes. 
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4 Retarded Green Functions and Criticality 



In this section, we first obtain formal expressions for the retarded Green functions of the vector current 
and energy momentum tensor operators of the boundary theory. By formal expressions we mean that 
we write formulas relating the above-mentioned retarded Green functions to IL±, to be introduced 
below. 

Consider first the solution of the equations (3.5), (3.6) and (3.7) asymptotically, for u — > oo. For 
each of the fields in this problem, there is a constant asymptotic value that we denote by a hat. The 
solutions take the form 



h\(u 
a y (u 



oo 



+ 



7T 



wh x +qh y f ) + 



w 

U 3 



+ ..., 



+ ..., 



OO 



- , % 
a y + — - 
y u 



3(c, 5 



vo- 



ii- 



(4.1) 
(4.2) 
(4.3) 



In these expressions, all quantities are understood to be functions of t» and q even if not explicitly 
indicated. Here, h x , h y t and d y are regarded, as usual, as the sources for the corresponding operators 
in the dual field theory, while Tr x y ,n y t and ir y are the corresponding one-point functions of those 
operators. The latter are understood to be functions of the sources, and that functional dependence, 
given the on-shell action, encodes the correlation functions of the dual field theory. 



The renormalized action in our case is given by 

6 



(Fx 



g[R + 



2K 



dM 



L 2 
4 

" L 



L 2 F^ 

d 3 x 

dM 



h\-L 



d 3 x 



dM 



|7| {3) ^, 



(4.4) 



where 7^„, K and are, respectively, the induced metric, the trace of the second fundamental 
form and the intrinsic curvature of the 3-dimensional boundary dM. The counterterm proportional to 
removes a linear divergence in the on-shell action for the Reissner-Nordstrom AdS4 background, 
with no contribution to its finite part. In order to extract Green functions that we are interested in 
holographically, it is enough to consider the (gravity) action up to quadratic terms in the fluctuations. 
One finds 



lim C 



J d 3 x [u A hy t hy' t - u A f{u)h x y h x ' y +Au 3 (l - fin)- 1 ' 2 ) (h\W> t - f(u)h x y h x y ) 



12 

u A f'(u)h x y h x y ^ (n a yh V t + u 2 f{u)a y a y ) 

A* 



- -^fiuy 1 ' 2 (d t h x y - d x hy t 



(4.5) 



where C = r /i 2 /( 12K 4)- 

The holographic prescription for evaluating the matrix of retarded Green functions when the op- 
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erators mix was given in [36] following [34, 35] 3 . Because we have rotational invariance in the spatial 
plane, we will always take the y-component of momentum to vanish, and hence it is most convenient to 
Fourier transform the x-direction only. Having imposed infalling boundary condition at the horizon, 
one writes the renormalized on-shell action in the form 

S .s. = J dy j ^| 0i(w,q)Jij(w,q)^(-w,-q), (4.6) 
where <pi = {h y t , h x y ,a y } in our case. The prescription of [36] then reads 

In momentum space, (4.5) evaluates to 

S ren = 3fi 2 C J dy J ^| [ 7 r :C j / (t»,q)^ J/ (-t»,-q) - 7^(10, -q) 

4 i 
+ -2 7r y (w,q)a w (-w,-q) + ... , (4.8) 

where the ellipsis contains terms which contribute at most to contact terms in two-point functions. 

Now, as we recalled above, the {ir X y,ir y t,Tr y } are regarded as functions of {h v t , h x y , a y }, and for 
present purposes, they can be regarded as linear functions. Also, by assumption, none of the cor- 
responding operators have one-point functions when the sources are removed. The vevs are not all 
independent in this case though. Indeed, the constraint (3.8) gives one more piece of information 

4 

C[lT X y + X»TT y t = -X»CL y (4.9) 

This equation is easily understood in the boundary theory. The diffeomorphism Ward identity [21] 
takes the form 

- 9ux(x)V»(f ""(z)> + F A(U (x)(>(x)) = 0. (4.10) 
whose y-component evaluates to (4.9) at linearized order given the source a y , and our notation. 

Thus, there are really only two independent vevs here, in accordance with the fact that the bulk 
problem can be reduced to the two gauge-invariant fields. What is more, the system is in fact diago- 
nalized by passing to the master fields. Asymptotically, we write the master fields as 

$ ± (u^oo)~$± ^1 + 5^ + ...^ , (4.11) 

where $-t and TT-|- are functions of to and q , and we have in mind that the master fields satsify infalling 
boundary conditions at the horizon. 



3 See also [37, 38] for a recent more systematic treatment of holographically computing Green functions in cases where 
the boundary theory operators mix. 
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By comparing asymptotic expansions, one easily deduces that 

$± = 3/xq (voh x y + q/^ t ) - 6 (l ± Vl+<l 2 ) ay, (4.12) 

while 

Mi = 3-71-% - 6 (l ± Vl+q 2 ) vr y , (4.13) 

where we have used (4.9) to simplify. It should be clear then that the on-shell action can be written 
as a function of the sources {h y t ,h x y ,a y } and II -t. The latter two functions (of no and q) must be 
determined numerically. The matrix of retarded two-point functions in the shear channel are then 
determined by these two functions. This is in fact consistent with the requirements of the Ward 
identities. 

Returning to the on-shell action, using the holographic prescription given above, we deduce the 
retarded correlators of J y , T xy and T y t 



G yt:yt (m,q) 



=/x 2 q 2 Gi(t»,q), 
= - // 2 qt»Gi(tt>,q) 
=^ 2 t» 2 Gi(t»,q), 



- -/xq 2 [Gi(tt>,q) + G 2 (tt>,q)], 



G xyjXy (tx>,<\) 
G yt , y (m,<\) 

G xy:y (m,q) =^/iqt» [Gi(w,q) + G 2 (n>,q)] , 
G y ,y(m,q) 



where 



Gi(t»,q) = 



G 2 (t»,q) = 



-4G 2 (t»,q), 

:!( r (i - A/T+q 2 ) n+ - (i + Vi + e) n- 
(1 + Vi + q 2 ) n+ - (i - TT+q 2 ) n_ 



3G 
V / l+q 



(4.14) 
(4.15) 
(4.16) 

(4.17) 

(4.18) 
(4.19) 

(4.20) 
(4.21) 



We have implicitly assumed here that the bulk equations of motion are to be solved by placing 
infalling boundary conditions on $±(u). It can easily be verified that the prescription (4.7) implies 
that G X y jyt (v>,<{) = Gy tjXy (vo,<\), Gy iXy (x0,c\) = G xy :V (t» , q ) , and up to a contact term Gy !yt {v>, q) = 
G yt , y (m,<\). 



4.1 Frequency Expansions of $ ± and II ± 

The expressions we derived above for the retarded Green functions are not very informative unless 
n_l_(t»,q) is known. One thing that can be done semi-analytically is to develop a series expansion in 
to. This is an extension of our earlier results to finite momentum. This analysis makes use of the fact 
that the near horizon geometry is AdS 2 x R 2 . 

Let us expand &±(u) in a power series in vo: 

*o±(«) = + w*8i(«) + «> 2 *S» + • ■ ■ , ( 4 - 22 ) 
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where the subscript denotes the outer region (see the discussion below). Since the black hole back- 
ground is extremal, taking the limit of to — > becomes somewhat more subtle near the horizon. This 
is due to the fact that in the extremal case f(u) has a double zero at the horizon. A way to handle 
this issue was explained in [17]. Basically, one realizes that near the horizon the equations for <3?-i- in 
(3.16) organize themselves as functions of ( = wq, where, using (2.9) and (3.4), we have 

u = l + ^L. (4.23) 

Since the coordinate £ is the suitable radial coordinate for the AdS2 part of the near horizon geometry, 
the tx> expansion of &± in that region can be written 

3j±(C) = *?2(C) + »*£(0 + » 2 *?2(C) + ■ ■ ■ , (4-24) 

where the subscript I denotes the inner region. Having implemented the (infalling) boundary condition 
at the horizon [34, 35] in the £ coordinate, one then matches the inner and outer expansions in the 
"matching region" where the £ — > and xo/C, — > limits are taken. Since the differential equations 
(3.16) are linear, and we also require that the solutions for the higher order terms in the expansions 
(4.22) and (4.24) do not include terms proportional to the zeroth-order solutions near the matching 
region, we just need to match 3>j±(C) to <E>q]_(u) near that region [17] at leading order. 

4.1.1 Inner Region 

Substituting (4.24) in (3.16), we find that the leading terms satisfy 

-*S"(0 + (2!±1±M+Z _ ij ,g (0 = o. (4 . 25) 

Equations (4.25) are identical to the equations of motion for massive scalar fields with effective AdS2 
masses 

™±L 2 2 = 1 + y ± \A+<t 2 - (4-26) 

Thus $^ (£) source operators 0± in the IR CFT with conformal dimensions <5± = v± + \ where 

v± = \ ^l + 4m 2 ± L 2 2 = X - ^5 + 2q 2 ± 4^1+^2. (4.27) 

Having imposed the infalling boundary conditions at ( — > oo, <3?^(e) take the following form near the 
matching region (where one takes the ( — > and vo/( — > limits) 

^fl(u -> 1) = \(u - l)--2 +u± + g±(w)(u - . (4.28) 

Note that in writing (4.28) we used (4.23) to express ( in terms of u. Also, we chose a specific 
normalization for Such a choice does not affect the calculation of the boundary theory Green 
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functions. In (4.28), G±(t») denote the retarded Green functions of the aforementioned IR CFT 
operators 0±, and are given by [17, 26] 



-2v±e 



-l1Tl/± 



r(l - V±) /VO\*v± 



(4.29) 



T(l + v±) \ 2 J 

The analytic structure of Q±{t») is the same as the retarded Green function of an IR CFT scalar 
operator discussed in [17]. Namely, for generic values of momentum q, there is a branch point at 
to = 0, and moreover the branch cut is chosen to extend along the negative imaginary axis in the 
complex tx>-plane. 



4.1.2 Outer Region and Matching 

In the outer region, the equations for $^]_(u) are obtained by setting t» = in (3.16). It is easy to 
show that near the matching region, the solutions for Qq±(u) are given by a linear combination of 
(u — l)~z +u± and (u — l)~^~ v± . For ease of notation, we define 



^ 0) (u) = («-l)-5+"± + ... ) 
Then, matching $jj±( u ) to (4.28), we obtain 

»(?) (,.\ - 



£i 0) («) = («-!)-§-* + ..., 



u 



1. 



(4.30) 



(4.31) 



The discussion here follows that of [17]. To higher orders, we can write 



rj±{u) = r]±\u) + vor] y ±>(u) + w' z r)^'(u) + . . . , 
£±(u) = ei 0) («) + to£\u) + tt) 2 ^) + . . . , 



2J2). 



(4.32) 
(4.33) 



where ??£ l>0 ^(u) and £± >0 \u) are obtained demanding that in the u — > 1 limit, they are distinct from 
rj±\u) and £±\u), respectively. Thus, 



$0±(«) = [r)±(u) + </±(n>)£t(u)] . 
Near u — > oo, one can expand r]±\u) and £±^(it) as follows (here n > 0) 





1 






li 


(n) 


1 




U 



(4.34) 

(4.35) 
(4.36) 



where the coefficients are all functions of q. So, asymptotically, we obtain 



a (°) + + o (to 2 ) 



do) 



+ £±(n>) ci 0) +tx>4 lj + O(t» 2 ) 



(1) 



&w+n>&w+0(» 2 ) + £t(») +»4 i; +0(w 2 ) 



(4.37) 
(4.38) 



By virtue of (4.14)-(4.19), these determine the frequency dependence of the retarded Green functions 
in the boundary theory, for small to. 
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4.2 Criticality: Emergent IR Scaling 



Although the Green functions (4.14)-(4.19) can be numerically computed by computing the coefficients 
a~± , b±^ , c*± , d±^ , it turns out that one can analyze the low frequency behavior of the Green functions 
without knowing these coefficients explicitly. For the cases of scalar and spinor operators, such analyses 
were performed in [17] where a variety of emergent IR phenomena were observed. In this subsection, 
we perform similar analyses of the retarded Green functions of the vector current and the energy- 
momentum tensor operators of the boundary theory at low frequency and determine their emergent 
IR behaviors. As we will see momentarily, one of the reasons which makes the low frequency analyses 
of the above-mentioned retarded Green functions non-trivial and interesting is the fact that there are 
now two IR CFT Green functions G±(b>) both of which can potentially play a role. Not surprisingly, 
the mixing of Q±(vo) in the expressions (4.14)-(4.19) for the boundary theory Green functions is due 
to the coupled nature of the electromagnetic and gravitational perturbations we are considering in the 
charged (extremal) black hole background. 

As defined in (4.27), v± are real implying that a±\ b± \ and d±^ are all real. Thus the complex 
part of Tl± comes entirely from the IR CFT Green functions Q±(yo). Suppose for generic momentum 
q, the product a^a^ is non- vanishing. Expanding (4.20) and (4.21) for small t», and noticing that 
2u + > 3 while 2z/_ > 1, one then deduces that Im Gi(t»,q) and Im G2(t»,q) exhibit scaling behavior 
at small frequency: 

ImGi(»,q)~ . (1 + Vl + q 2 ) e (q) Im Q- (») [1 + . . .] oc m 2u - , (4.39) 

V 1 + < » v 
3C 



Im G 2 (w, q) ^ , ( 1 - 7T+q2 ) e (q) Im Q- (») [1 + . . .] oc xa 2v - , (4.40) 

V 1 + < » 



where 



jio) / c (o) d m\ 

Substituting (4.39) and (4.40) into (4.14)— (4.19), the spectral functions (the imaginary part of the 
aforementioned retarded Green functions) of the vector current J y and the energy-momentum tensor 
operators T y t and T xy of the boundary field theory show scaling behavior at small frequency 

Im G y tyt(to,n) oc vo 2u -, ImG I „(ro,q)ocn) 2+2,/ -, Im G ytxy (vo, q) oc vo l+2u - , 

Im G y ty(t»,q) oc to 2 "-, Im G xyy (vo, q) oc vo l+2v - , Im G yy (t», q) oc rx> 2u - . 

Note that the low-frequency scaling behavior of the spectral functions above is emergent meaning 
that it is a consequence of the fact that the near horizon region of the background geometry (which 
translates into the IR physics of the boundary theory) contains an AdS 2 factor. Although eo depends 
on details of the background geometry, the scaling behavior vo 2v - does not change by changing the 
geometry in the outer region as long as the AdS 2 part of the near horizon geometry is kept intact. 
Also, note that because v± are strictly real, the retarded Green functions (4.14)-(4.19) do not exhibit 
the log-periodicity behavior observed in [16, 17] for the cases of charged scalar and spinor operators. 

Of course, it is also of interest to consider the physical poles and branch points of the correlation 
functions at arbitrary q and rx>. We turn to a study of that problem in the following sections. 
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5 Retarded Green Functions and Quasinormal Modes 



The expressions (4.14)-(4.19) for the retarded Green functions of the boundary field theory indicate 
that they generically have poles whenever <&±(t»,q) = 0. We denote by t»^(cj) the frequencies at which 
*5±(t»,C|) vanish. In order to obtain our boundary theory Green functions holographically, we imposed 
infalling boundary condition on &± at the horizon. By further imposing &± to vanish asymptotically 
we are essentially solving for their quasinormal modes (QNMs) in the extremal Reissner-Nordstrom 
AdS4 black hole background. 

In the context of the AdS/CFT correspondence, the connection between the QNMs of gravitational 
backgrounds and the singularities of the retarded Green functions of the dual boundary theories was 
first suggested in [39]. It was then argued in [34] that the quasinormal frequencies of the decoupled 
fluctuations of a black hole background coincide with the poles in the retarded Green functions of 
the corresponding dual operators in the boundary theory. The authors of [40] then argued that 
in cases where the fluctuations are coupled, in order for the quasinormal frequencies to coincide 
with the poles in the retarded Green functions of the dual operators, one should consider the gauge 
invariant combinations as the right set of variables and impose on them infalling boundary condition 
at the horizon together with an asymptotic normalizable boundary condition. In our case, as it is 
seen from the linearized Einstein-Maxwell equations (3.5)-(3.8), the electromagnetic and gravitational 
perturbations of the background are certainly coupled. In section 3 we constructed the gauge invariant 
combinations X(u) and Y(u) in terms of which the linearized Einstein-Maxwell equations are still 
coupled. We decoupled the equations for the gauge-invariant modes by introducing new variables 
&±(u). These variables are themselves gauge invariant by virtue of their definition (3.15) in terms 
of X(u) and Y(u). Thus, with $±(it) being gauge invariant and satisfying decoupled second-order 
differential equations, we can analyze their QNMs. 

To our understanding of the literature, the electromagnetic and gravitational QNMs of the extremal 
Reissner-Nordstrom AdS4 black hole have not been studied before 4 , which makes finding these QNMs 
an interesting problem in its own. Ref. [45] is a nice review of QNMs from both general relativity and 
the AdS/CFT perspectives. 

There are only very few backgrounds whose QNMs are known analytically. Although there are 
semi-analytic methods for calculating small or large (magnitudes of) quasinormal frequencies, one is 
usually forced to do numerics in order to extract the generic values of such frequencies. We refer 
the reader to [45] for an extensive list of analytic and numerical methods available in the literature 
on how to compute the QNMs of different backgrounds. In this section, we numerically compute 
the quasinormal frequencies of the (shear channel) electromagnetic and gravitational perturbations of 
the extremal Reissner-Nordstrom AdS4 black hole. Along the way, we find it useful to compare our 
zero-temperature numerical results with the analogous QNMs of the non-extremal Reissner-Nordstrom 
AdS4 black hole background. Although, the results in the non-extremal case are not new, we present 
them here for the purpose of comparison and to study the T — > limit (at fixed fi). By doing so, we 

4 The QNMs of the four-dimensional asymptotically flat extremal Reissner-Nordstrom black hole were studied in 
[42]. Ref. [43] analyzed the scalar, electromagnetic and gravitational QNMs of non-extremal Reissner-Nordstrom AdS4 
black hole. Recently, the authors of [44] studied the QNMs of a massive, charged scalar field in the background of the 
electrically, as well as the dyonic extremal Reissner-Nordstrom AdS4 black hole. 
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will see more clearly how the analytic structure at T = is related to that of finite temperature. 



5.1 Generalities 

Before numerically computing the quasinormal frequencies for our problem, we discuss some general 
features that one expects to observe. 

5.1.1 Branch Cut 

For small to, the analytic structure is determined by the small frequency analysis that we presented 
earlier, pertaining to the near-horizon geometry. For generic values of q, Q±(w) are multi- valued and 
the boundary theory Green functions are then generically multi-valued as well. When q = 0, 2u± are 
integers: indeed 2v+ = 3 and 2z^_ = 1. In this case, the multi- valuedness of the Green functions is 
due to the existence of terms like t» m logt» in their expressions, where m is a positive integer, in the 
extremal case. 5 For example, one can show that at small vo 

ivo 

G xv q = 0) oc , (5.1) 

y ' yy ' ; 1 + ci» log w + c 2 itx> y 1 

where c\ and C2 are numerical factors. Thus, the QNMs of <&± will also be multi-valued (in the 
extremal case). There is another, perhaps, more direct argument for the existence of the branch cut in 
the QNMs of 3>±, which is based on analyzing the "tail" of the 3>± fluctuations in the extremal black 
hole background [46] . An argument of this type was used in [44] to show the existence of the branch 
cut in the QNMs of a chargeless massless scalar in the extremal Reissner-Nordstrom AdS4 black hole 
background, and can easily be adopted to our case, as well. To that end, define the (dimensionless) 
tortoise coordinate u* via 

7T = -^TY < 5 ' 2 > 
au u z j(u) 

The asymptotic boundary of the background in the tortoise coordinate is obtained by taking u* — > 
whereas u* — > — oo is the near horizon region. Using the tortoise coordinate, we can write the equations 
in (3.16) in the form of Schrodinger equations (where we think of q as being fixed) 

$±(u*) + [3o> 2 - V±(u*)] $±(u*) = 0, (5.3) 

where the potentials V±(u*) are given by 

V±M = (l ± VT+^) + 3q 2 /W - u 3 f(u)f(u). (5.4) 

It is understood that in the above expressions for V±(u*), one should use (5.2) to express u in terms 
of u*. Near the horizon, one easily finds 

V±(ii* -)• -oo) ~ (5.5) 



5 As we will recall below, such logs are absent at finite temperature, and the cut resolves into a series of Matsubara 
poles along the negative imaginary axis. 
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Such a power law behavior in the asymptotic form of the potentials is a characteristic of the background 
being extremal (as opposed to a non-extremal background for which the potential dies off exponentially 
as it* —7- — oo). It is this power law behavior which gives rise to a branch cut [46] (for small frequencies), 
extended in the negative imaginary axis, for the QNMs of &±. We will shortly see the appearance of 
this branch cut in our numerical plots of the QNMs. 



5.1.2 Stability 

Suppose <I>-|- are QNMs with frequencies n>^. So, Q± are ingoing at the horizon and vanish asymp- 
totically. Following [47], we will see shortly that, using the properties of the potentials V±, these two 
boundary conditions constraint t»^. Since &± are infalling at the horizon, one has e~ ty ^ WnT $± ~ 
e -i v / 3t» n (r+u») w h ere r _ rot/L 2 is the dimensionless time coordinate and u* is the tortoise coordinate 
defined in (5.2). Note that r + u* is the dimensionless ingoing Eddington coordinate. Define 

$±(u„) = e" i>/Sw » u 'tf±(u,). (5.6) 
Substituting (5.6) into (5.3), and writing the result in terms of the coordinate u, we obtain 

[« 2 /(u)*±(«)]' - 2iV3w^'±(«) - U±(u)V±(u) = 0, (5.7) 
where the potentials U±(u) are given by 

U±(u) = A (i ± TTT^) + ^q 2 - uf(u). (5.8) 



We multiply (5.7) by ^±(u) and integrate the result to obtain 

r°° r 

[u 2 f(u)y±(u)^' ± (u)}™ - / du u 2 f(u)\^' ± (u)\ 2 + 2iV3m±^±(u)&±(u) + U±(u)\y±(u)f 



= 
(5.9) 

where we have also performed an integration by parts. Using (5.6) and (4.11) and the fact that <&±(u) 
are assumed to be QNMs, one obtains that as u — > oo, u 2 ^'±(u) is finite and ^±(u) vanishes. At the 
horizon, both *f!±(u) and ^±(u) are finite and f(u) vanishes. Thus, all together, one concludes that 
[u 2 f (u)^ ±(u)^' ± (u)]'^ ) = 0. Take the complex conjugate of (5.9) and subtract the result from (5.9) 
to obtain 

/oo 
du [w±*±(«)*' ± («) + w±*±(u)*' ± («)] = 0. (5.10) 

Integrating by parts the second term in (5.10), and noting that i &±(u) vanishes at u = oo, we obtain 

/oo • 
d«*±(«)*' ± («) = |tf±(u = 1)| 2 . (5.11) 

Since ^>±(u = 1) is finite, if Re vo^ / equation (5.11) then implies that Im vo^ / 0. Substituting 
(5.11) into (5.9), we obtain 



L du 



u 2 



f(u)\*' ± (u)\ 2 + U ± (u)\y±(u)\ 2 +v / 3(lmtx>±) |w± | 2 |tf±(« = 1)| 2 = 0. (5.12) 
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For all (real) values of q, U+(u) is always positive definite. This is because in order for U+(u) to be 
always positive definite, one has to have 3q 2 + 2^1 + q 2 — 2 > 0, which is satisfied for all (real) values 
of q. Thus, the integral (with subscript +) in (5.12) never becomes negative. From equation (5.12), 
we then reach the conclusion that Im t»+ < 0. This result indicates that the <J?+(u) fluctuation does 
not cause instability, as e~ t ^ 1x>nT & + decays in time for all u> + . For U-(u), on the other hand, it can 
easily be shown that the condition for it being positive definite turns out to be 3q 2 — 2-y/l + q 2 — 2 > 
which is true only if q 2 > Therefore, for |q| > |, the integral (with subscript — ) in (5.12) is 
positive definite yielding Im t»~(q) < 0. For |q| < |, there are ranges of u for which U-{u) becomes 
negative. As a result, when |q| < |, without knowing the solution equation (5.12) cannot 

by itself determine the sign of Im t»~(q). However, as we will see later on, it is easy to numerically 
confirm that even for |q| < | one obtains Im t»~(q) < 0, so the <£_(«) fluctuation does not cause 
instability either. 

In the background geometry we are considering, the quasinormal frequencies of <E>± are further 
constrained. Using equations (5.3) together with the infalling boundary condition (5.6), one easily 
concludes that if t»^(q) is a quasinormal frequency of Q±, so is — t»^(q), that is to say that the 
quasinormal frequencies appear as pairs with the same imaginary part but opposite real part. On the 
field theory side, the pairing of the poles of the retarded Green functions (4.14)^(4.19) is due to the 
fact that the boundary theory is parity invariant. 6 



5.2 Numerical Analysis: Matrix Method 

To compute the quasinormal frequencies of <&± (in the extremal case) numerically, we use the "matrix 
method" of [48]. Unlike the algorithm of [47] which is known not to work for extremal black hole 
backgrounds 7 , the matrix method of [48] is believed to be applicable here. The authors of [44] used 
this matrix method to find the quasinormal frequencies of a charged minimally-coupled scalar field in 
the extremal Reissner-Nordstrom AdS4 background. In what follows we first present our numerical 
results for the quasinormal frequencies of $± in the extremal case, then compare their general features 
with the analogous quasinormal frequencies when the Reissner-Nordstrom AdS4 black hole background 
is non-extremal. 



5.2.1 Zero Temperature 

For numerical purposes it is more convenient to switch to a new radial coordinate z = 1/u. In this 
coordinate, z = 1 represents the horizon while z = is the asymptotic boundary. We write the 
equations (3.16) in terms of the new coordinate z and arrive at 



[f(z)*'±(z)]' + \-zf'(z) + |^ (tx) 2 - /(z)q 2 ) - 6z 3 (l ± yi+q 2 ) 



$±(s) = 0. 



(5.13) 



To apply the matrix method of [48] , we first isolate the leading behaviors of <3?-t (z) at the horizon and 
the boundary. The rest is then approximated by a power series in z around a point zq such that the 

6 This symmetry of the quasinormal modes would be broken, for example, in the presence of a magnetic field. 
7 This is because for extremal black holes, the horizon is an essential singularity. 
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radius of convergence of the series covers both the horizon and the boundary. That point in our case 
is taken to be zq = \. Thus, we write 



m=0 ^ ' 



(5.14) 



Substituting (5.14) into (5.13), we obtain a set of M+l linear equations for M+l unknowns a p (t»,q)'s, 
and write them in a matrix form 



M 



Yj A m P (o>,q)aJ(t»,q) = 0. 



(5.15) 



m=0 



Suppose now we choose a specific value for the spatial momentum q. The quasinormal frequencies vo^ 
for that particular value of q are solutions to the following equation 

det A± p (w±, q) = 0. (5.16) 

Since <I>-|- in (5.14) are approximated by power series, the more terms kept in the series, the more 
accurate results one obtains for the quasinormal frequencies of &±. 



(a) 



(b) 




a. 
3 




Figure 1: Electromagnetic and gravitational quasinormal frequencies (in the shear channel) of extremal Reissner- 
Nordstrom AdS4 black hole. Plot (a) shows the quasinormal frequencies of <J>+ and plot (b) shows the quasinormal 
frequencies of For both plots, q = 1 and M — 100. As M is increased the poles on the negative imaginary axis get 
closer to one another and form a branch cut in the M — > 00 limit. 



The plots (a) and (b) in Figure 1 show respectively the quasinormal frequencies of <E> + and 3>_ for 
q = 1, where, due to space limitations, we have only shown a handful of them. We used M = 100 
for both plots. It can easily be verified numerically that, for both plots, as M is increased the poles 
on the negative imaginary axis become closer to one another. This observation suggests that the 
sequence of discrete poles on the negative imaginary axis is due to approximating $-t in (5.14) by 
power series. In particular, although not numerically accessible, it should be the case that in the 
M — > 00 limit there is actually a branch cut [44]. This picture is consistent with the two arguments 
we presented earlier for the existence of the branch point at the origin in the extremal case. Notice 
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that the plots are symmetric under xo^ — > — as they should be. One also sees that the off-axis 
poles are approximately equally spaced and line up on almost straight lines. Aside from the poles on 
the negative imaginary axis, which coalesce as M is increased, the qualitative behavior of the plots 
does not change by increasing M. In particular, we do not see any poles crossing into the real axis, 
indicating that these QNMs do not cause instabilities in the system. This observation is consistent 
with the argument we gave earlier in subsection 5.1.2. The argument there was not powerful enough 
to determine the sign of Im t»~(q) for q £ (— |, |). The plot (b) of Figure 1 (for which q = 1) shows 
that the imaginary part of v»~ is negative. One can also plot the quasinormal frequencies of <3>_ for 
other values of q in the range (— |, |) and observe that Im xx>~ is always negative. Indeed we have 
shown one such plot in Figure 2(b). 

The plots in Figure 2 show the quasinormal frequencies of &± for q = 1/ y/3. Although not shown, 
we considered other values of q and were able to verify that the qualitative patterns of the plots shown 
in Figures 1 and 2 stay the same as q is varied. Excluding the poles on the negative imaginary axis, 
which form a branch cut in the M — > oo limit, we observed that independent of what value q takes, 
there are no small quasinormal frequencies for <&_, while there are only two for $ + (which are the 
ones closest to the real axis). Here, by small quasinormal frequencies we mean those which satisfy 
|t»^| < 1. Therefore, one can conclude that at small frequencies |t»| < 1, there are only two single 
poles in the boundary theory retarded Green functions (4.14)-(4.19). Moreover, those poles do not 
approach zero as q — > 0. 



(a) 



(b) 



a. 




a. 

"J 




Figure 2: Quasinormal frequencies of <2>± (in the extremal case) for q = 1/V3- (a) Quasinormal frequencies of 3>+. (b) 
Quasinormal frequencies of M = 100 for both plots. 

The analysis at q = is quite a bit simpler, but is qualitatively similar. In this case, the relevant 
equations are 



z 2 f(z)h*'>(z) + [z 2 f'(z) - 2zf(z)] h%(z) + j^z 2 h\(z) = 0, 
f(z)a'i(z) + f{z)a' y {z) + - 12z 2 ) a y (z) = 0, 



(5.17) 
(5.18) 



which are obtained from equations (3.13) and (3.14) by changing the radial coordinate to z = 1/u. 
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(Note that h x y {z) is gauge invariant at q = 0, while a y (z) is gauge invariant for all q.) We can find 
the quasinormal frequencies of h x y {z) and a y {z) similar to the steps we performed above to obtain the 
quasinormal frequencies of <&±. The only difference is that the normalizable mode of h x y {z) falls off 
asymptotically as z 3 . Therefore in writing the analog of (5.14) for h x y {z) the factor of z should be 
replaced by z 3 . The results are shown in Figures 3(a) and 3(b). Note that equation (5.17) is the same 
as the equation of motion for a minimally-coupled chargeless massless scalar in the background of the 
extremal Reissner-Nordstrom AdS4 black hole. The quasinormal frequencies of this scalar have been 
numerically computed in [44] using the same matrix method and a (density) plot resembling the one 
in Figure 3 (b) was generated there. 

(a) (b) 




-4 -2 2 4 -4 -2 2 4 

Re (ui/iX) Re (ai/pi) 

Figure 3: Quasinormal frequencies for q = of a y and h x y in the extremal Reissner-Nordstrom AdS4 black hole 
background. Plot (a) shows the quasinormal frequencies of a y and plot (b) are the quasinormal frequencies of h x y . We 
set M = 100 for both plots. 

The pattern we observe in the two plots of Figure 3 is very much similar to the ones in Figures 1 
and 2. The plots were obtained for M = 100, but the pattern stays almost the same for higher values 
of M. One still observes that the poles on the negative imaginary axis coalesce while the locations of 
the poles on each side of the imaginary axis do not change much. Also, we neither find poles crossing 
into the upper half plane, nor do we find poles which vanish. As it is seen from Figure 3(b), there are 
no (diagonally-oriented) quasinormal frequencies for h x y at small \rx>\ < 1 frequencies while there are 
two such frequencies for a y in that range. The quasinormal frequencies of a y in Figure 3(a) can be 
used to estimate, along the lines of [44], the UV behavior of a one-loop contribution to the free energy 
of the extremal background by the gauge field a y . Knowing the one- loop corrections to the free energy 
enables one to compute a "1/jV correction" to charge conductivity. 

5.2.2 Finite Temperature 

It is instructive to compare our results for the quasinormal frequencies of §± in the extremal case to 
the quasinormal frequencies of 3>± in the background of non-extremal Reissner-Nordstrom AdS4 black 
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hole. Written in terms of the z-coordinate, equations (A. 9) take the form 



z 2 [/(*)*' ± (*)]' + 



*±(*) = o. 



(5.19) 



where Q 2 < 3 in the non-extremal case. The definition of &± at finite temperature is given in (A. 7). 

(a) (6) 





Figure 4: Quasinormal frequencies of $± in non-extremal Reissner-Nordstrom A0IS4 black hole background. Plot (a) 
shows the quasinormal frequencies of "I>+ while plot (b) shows the quasinormal frequencies of For both plots, q = 1, 
M = 100, and T = 0.09/it. 



In the non-extremal case, we also use t»„ to denote the quasinormal frequencies of Q±. In order to 
find rx>^ we use the matrix method of [48] and follow the exact same steps as we did in the extremal 
case except that (5.14) is now replaced by 



M / 1\ 

M*)=/(zr ro/r *E°m( ro >*)(*-2 

m=0 ^ ' 



(5.20) 



where we have defined r = AirT/ n, with T and fi being the temperature and chemical potential, 
respectively. The finite-temperature quasinormal frequencies of &± have been plotted in Figure 4. 
Only a handful of n>^ have been shown, though. To generate the plots, we chose a temperature of 
T = 0.09/i. Also, the plots were obtained for q = 1 and M = 100. The branch cut in Figures 1(a) 
and 1(b) are replaced at finite temperature by a sequence of almost equally-distanced poles along 
the negative imaginary axis. Note that the poles are all in the lower half of the frequency plane, 
indicating that they do not give rise to instabilities. As q — > 0, none of the quasinormal frequencies 
of <!>+ approach zero. On the other hand, the lowest quasinormal mode of shown in Figure 4(b) 
by a larger red dot, approaches zero as q — > 0. This is the leading hydrodynamic pole. Thus, one 
concludes that the finite temperature generalization of (4.14)-(4.19) all have the same hydrodynamic 
pole, namely the diffusion constant T> is the same for all of them. In Figure 5, we plotted the lowest 
quasinormal frequency of <&± (shown in Figure 4(b) by a larger red dot) as a function of q, and fitted 
the plot by a quadratic function. The solid black curve on the plot shows the fit. From the slope of the 
fitted curve, we obtain V ~ 0.17///, with T = 0.09/i. Substituting the numerically obtained diffusion 
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Figure 5: The lowest quasinormal frequency of <E>_ as function of q. The temperature is T = 0.09/1, and M = 100. The 
solid black curve is the quadratic fit: ivo ~ 0.17q 2 . 



constant into the hydrodynamic equation r/ = (Ts + /Up)X>, where p, s and p are respectively given in 
(2.4) and (2.6), one easily obtains rj/s ~ 1/Att. 8 



(a) 



(b) 





Figure 6: (a) Finite-temperature quasinormal frequencies of a y at q = 0. (b) Finite-temperature quasinormal frequen- 
cies of h x y . The temperature is T = 0.09/x, and M = 100. 



We present the case of q = separately. The relevant equations in this case read 
z 2 f(z)h x ;(z) + [z 2 f(z) - 2zf(z)} h x ' y {z) + Cf—**h x v (z) = 0, 

4z 2 ) a y (z) = 0, 



f(z)al(z) + f(z)a' y (z) + Q 2 ^ 



(5.21) 
(5.22) 



which are equations (A. 5) and (A. 6) expressed in terms of the coordinate z = 1/u. We find the finite- 
temperature quasinormal frequencies of h x y {z) and a y {z) following what we did above to obtain the 



8 The sign ~ is because of numerical uncertainties. One can of course show analytically that the ratio should be equal 
to 1/4tt [28, 29]. 
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finite temperature quasinormal frequencies of &±. The only difference is that in writing the analog 
of (5.20) for h x y (z) we replace the factor of z by z 3 . This is because asymptotically the normalizable 
mode of h x y (z) falls off as z 3 . 

The quasinormal frequencies of a y and h x y have been plotted in Figure 6 for a temperature of 
T = 0.09/x using M = 100. We see that, as for q / 0, there is no branch cut at finite temperature: it 
is replaced by a sequence of almost equally-distanced poles along the negative imaginary axis. Also, 
as it is seen from the plots, the finite-temperature quasinormal frequencies of a y and h x y stay away 
from the real axis. To the extent that we have checked numerically, this behavior of the quasinormal 
frequencies is unchanged as M is increased. 

There does not appear to be a qualitative difference in the quasinormal spectrum between the 
T = results and the finite temperature spectrum for small temperatures. To see that more clearly, 
we consider the limit T — > 0, holding fi fixed. For numerical simplicity, we have done this at q = 0, but 
do not expect an appreciable difference at finite q. Referring to Figure 6, we consider the following 
three quantities. We denote by d\ the difference between the smallest quasinormal frequency on the 
negative imaginary axis and the origin, by efe the distance between the first two poles on the imaginary 
axis and by the distance between the two off- axis poles in the first quadrant closest to the origin. In 
Figure 7(a), we demonstrate that as T is lowered, holding fi fixed, both d\ and di approach zero, while 
d,3 remains finite. Approaching T — > requires ever-more numerical precision, but the plots suggest 
that d\ and di extrapolate to zero at T = 0. This is consistent with the on-axis poles becoming the 
branch cut seen in the T = theory, as one might expect. 

(a) (b) 




0.00 0.05 0.10 0.15 0.20 0.00 0.05 0.10 0.15 0.20 



Tin Tin 

Figure 7: (a) Plots of d\ (red), d,2 (blue) and dz (green) as function of T//x for a y . (b) Plots of di, (fe and ^3 as 
function of T/fi for h x y . In the range of T/(j, £ [0.017,0.091], the red and blue curves almost overlap. The plots have 
been generated with M = 150. 

6 Conclusions 

In this paper, we considered a boundary field theory dual to the extremal Reissner-Nordstrom AdS4 
black hole, and analyzed the retarded two-point functions of the charge vector current and energy- 
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momentum (tensor) operators at finite frequency and momentum. The operators J y , T yt and T xy whose 
two-point functions we analyzed were in the shear channel of the boundary field theory, whose dual 
bulk modes a y , h y t and h x y were odd under parity, y — > — y. We found that for generic momentum, the 
retarded correlators of the aforementioned vector and tensor operators exhibit emergent scaling behav- 
ior at low frequency. The operators J y , T yt and T xy gave rise in the IR to two sets of scalar operators 
0± with conformal dimensions 5±. The low frequency scaling exponents of the spectral densities in 
the boundary theory depended only on <5_ (or equivalently, z/_). We explored the analytic structure of 
the retarded correlators for finite momentum and small frequency, and argued that there should exist 
a branch cut in these correlators at the origin, as well as a series of metastable modes corresponding 
to isolated poles in the lower half of the complex frequency plane. Taking linear combinations of the 
dual bulk modes, we constructed the gauge invariant modes from which we were able to numerically 
compute the spectrum of the (shear-type) electromagnetic and gravitational quasinormal frequencies 
of the background. The numerical computations were in agreement with what we argued, analytically, 
for the existence of the branch cut in the retarded correlators as well as their poles being in the lower 
half of the complex frequency plane. Turning on temperature, we numerically observed how the branch 
cut dissolves into a series of poles on the negative imaginary axis, and determined the corresponding 
dispersion constant of the leading pole, in agreement with previous results. We performed the shear 
channel analyses in the (2 + l)-dimensional boundary theory. One can straightforwardly generalize the 
same type of analyses for higher dimensional boundary theories dual to extremal Reissner- Nordstrom 
AdS(i+i black hole, and, not surprisingly, reach the same conclusions as we did for the case of the 
(2 + l)-dimensional boundary theory. 

Although more involved, analyzing the retarded correlators of the operators in the sound channel 
of the (2 + l)-dimensional boundary theory can similarly be performed. In the bulk, one deals with the 
modes which are even under the parity, y — >■ —y: ha, ht x , h xx , h yy , at, a x . Taking linear combinations 
of these modes, one can construct two gauge-invariant modes, and compute, either numerically or 
analytically for small t» and q , the spectrum of their quasinormal frequencies (having diagonalized the 
system of equations) . The sound wave dispersion relation t» = to (q ) can then be determined from the 
poles of, say, Grtt,tt(°>>q) at small tx> and q. It is interesting to verify that t»(q) — > as q — > 0. The 
analysis for the sound channel of the (2 + l)-dimensional boundary will be presented elsewhere. 
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A Finite Temperature Equations 



In the background of non-extremal Reissner-Nordstrom AdS4 black hole, the linearized Einstein- 
Maxwell equations for the (odd-parity) modes h y t , h x y , a y take the form 



f(r) r 4 /^'/(r) + 4r 3 /i^(r) + AQr 2 L 2 a' y {r) - ukL 4 h x y (r) - k 2 L 4 h y t (r) = 0, 
f(r) \r A f{r)h x l{r) + [r 4 f'(r) + 4r 3 f(r)] h x ' y (r)] + u 2 L 4 h x y {r) + ukL 4 h y t {r) = 0, 
f(r) [r 4 f(r)L 2 a;(r) + r 2 L 2 [r 2 f'(r) + 2r/(r)] a' y (r) + Qr 2 r 2 h?\(r) 
+ L 6 [u?-j(r)h 2 \ a y (r) = 0, 
Also, the constraint equation (the yr-component of the linearized Einstein equations) reads 

- 4 - • h y ' t {r) + r A f(r)k h x ' y {r) + AQr 2 L 2 oo a y (r) = 0. 



r co 



(A.l) 
(A.2) 

(A.3) 
(A.4) 



In the above expressions, f(r) is given in (2.4) and has a single zero at the (outer) horizon r^. Note 
that M = 1 + Q 2 . 

In the non-extremal case, we choose to rescale oj and k by ji. The dimensionless frequency and 
momentum, to and q, are defined in (3.4). Also, like the extremal case, we work with the dimensionless 
radial coordinate u = r/r$. The gauge- invariant combinations are given by the same expressions as in 
(3.11) and (3.12): X(u) = <\h y t {u) + voh x y (u) and Y(u) = a y (u). 

At q = 0, h x y {u) becomes gauge- invariant. Also, the linearized Einstein-Maxwell equations (A.l)- 
(A.4) give, in this case, two decoupled equations for h x y {u) and a y (u) 



u 4 f{u)h x l{u) + [u 4 f'(u) + 4u 3 f(u)] h x ' y {u) + Q 2 ^-h x y (u) = 0, 

J\ u ) 



u 2 /(n)a» + [u 2 f{u) + 2uf(u)] a' y {u) + 



Q 2 \ w 2 



u 2 f(u) u 2 



a y {u) = 0. 



(A.5) 
(A.6) 



For finite (non-zero) q, on the other hand, it can easily be worked out that the linearized Einstein- 
Maxwell equations yield two coupled second-order differential equations for X(u) and Y(u). Similar 
to the extremal case studied in the bulk of the paper, these equations can be decoupled. We first 
define 

2f(uW 



<f>±{u) 



<\f(n)u 3 , 2Q 2 



rx> 2 — /(u)q 2 " 



u [no 2 — f(u)q 2 



+ ug±{q) 



Y(u), 



where 



9±(q) 

The decoupled equations then read 
[u 2 f(u)& ± (u)]' 




1± 4 1 + 



16 



1 + 



Q 2 



-2 



(A.7) 



(A.8) 



uf'(u) + (n> 2 - /(u)q 2 ) - (q ) 



u 2 f(u) 



u° 



$±(u) = 0. 



(A.9) 
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